NASA TN D-5863 

■ e / 


NASA TECHNICAL NOTE 


CO 

SO 

oo 


a 




AFWL (WLOL) 
KIRTTi-AND AFB^ N MJ 


<c 


TURBULENT-BOUNDARY-LAYER 
HEAT-TRANSFER AND TRANSITION 
MEASUREMENTS WITH SURFACE 
COOLING AT MACH 6 

by Aubrey M» Cary, Jr. 

Langley Research Center 
Hampton, Va. 23365 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. • JUNE 1970 


TECH UBRARY KAFB, NM 




TECH LIBRARY KAFB, NM 





1. Report No. 2. Government Accession No. 

NASA TN D-5863 

3. Recipient's Catalog No. 

4. Title and Subtitle 

TURBULENT-BOUNDARY-LAYER HEAT -TRANSFER AND 
TRANSITION MEASUREMENTS WITH SURFACE COOLING 
AT MACH 6 

7. Author(s) 

Aubrey M. Cary, Jr. 

5. Report Date 

June 1970 

6. Performing Organization Code 

8. Performing Organization Report No. 

L-6825 

10. Work Unit No. 

129-01-20-09 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23365 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15^ Supplementary Notes 


The material presented herein is based on a thesis entitled "Turbulent Boundary Layer 
Heat Transfer and Transition Measurements With Extreme Surface Cooling at Mach 6" sub- 
mitted in partial fulfillment of the requirements for the degree of Master of Aerospace 
Engineering, University of Virginia, Charlottesville, Virginia, August 1969. 

16. Abstract 

This experimental investigation shows the effect of wall cooling on turbulent heat- 
transfer and boundary-layer transition for a sharp -leading- edge flat plate at Mach 6 and 
Reynolds numbers as high as in the Langley 20-inch hypersonic wind tunnel. Decreasing 
the ratio of wall temperature to total temperature ^0.19 ^ T^/T^ = 0.7^ had little effect on the 
turbulent heat-transfer coefficient. The Spalding and Chi method best predicted both the level 
and trend of the heat -transfer data with wall cooling for the virtual origin located at peak 
heating and using the modified Karm^ Reynolds analogy. The transition Reynolds number 
increased (30 to 60 percent) when the ratio of wall temperature to total temperature decreased, 
but no transition reversal occurred. 


17. Key Words (Suggested by Author(s) ) 


18. Distribution Statement 



Turbulent boundary layer 
Turbulent heat transfer 
Boundary-layer transition 
Hypersonic 
Surface cooling 


Unclassified 

- Unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

48 

$3.00 


*For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 



I 


TURBULENT-BOUNDARY-LAYER HEAT-TRANSFER AND 
TRANSITION MEASUREMENTS WITH SURFACE 
COOLING AT MACH 6^ 

By Aubrey M. Cary, Jr. 

Langley Research Center 

SUMMARY 

This experimental investigation shows the effect of wall cooling on turbulent heat- 
transfer and boundary-layer transition for a sharp-leading-edge flat plate at Mach 6 and 
Reynolds numbers as high as in the Langley 20-inch hypersonic wind tunnel. 

Decreasing the ratio of wall temperature to total temperature ^0.19 S T-^y/T^ i 0.?j had 
little effect on the turbulent heat -transfer coefficient. The Spalding and Chi method best 
predicted both the level and trend of the heat -transfer data with wall cooling for the virtual 
origin located at peak heating and using the modified Karman Reynolds analogy. The tran- 
sition Reynolds number increased (30 to 60 percent) when the ratio of wall temperature to 
total temperature decreased, but no transition reversal occurred. 

INTRODUCTION 

Turbulent boundary layers have been observed over large areas of flight vehicles 
up to Mach 6 (ref. 1) and may be expected at even higher Mach numbers on proposed large 
vehicles. Since the structural and thermodynamic design of such vehicles will depend 
upon knowledge of the surface heating and friction, prediction methods that allow an accu- 
rate description of laminar, transitional, and turbulent boundary-layer flow are required. 

Because of the difficulties inherent in any theoretical approach to calculate the turbu- 
lent boundary layer, numerous semiempirical methods for predicting turbulent boundary- 
layer characteristics have evolved. In several recent attempts (for example, refs. 2, 3, 
and 4) to establish the validity of available theories and methods for turbulent flow by com- 
paring predictions with experimental skin-friction and heat -transfer data, certain semi- 
empirical methods were found to be superior. Semiempirical methods most commonly 
found in the literature are the reference temperature (or enthalpy) methods (refs. 5 to 7), 

^The material presented herein is based on a thesis entitled "Turbulent Boundary 
Layer Heat Transfer and Transition Measurements With Extreme Surface Cooling at 
Mach 6" submitted in partial fulfillment of the requirements for the degree of Master of 
Aerospace Engineering, University of Virginia, Charlottesville, Virginia, August 1969. 
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the Winkler and Cha method (ref. 8), and the recent methods of references 2 and 9. Avail- 
able experimental data guided the formulation of each of these methods. 

Since most hypersonic flight vehicles will have low wall temperatures in comparison 
with the total temperature of the flow, prediction methods for turbulent flow must be valid 
in the low range of the ratio of wall temperature to total temperature ^Tw/Tt). Studies by 
Spalding and Chi (ref. 2) and Bertram and Neal (ref. 4) have shown a paucity of experi- 
mental data for flat-plate turbulent flow with ratios of wall temperature to total tempera- 
ture below 0.5. In addition, reference 4 showed that predictions of h 3 ^ersonic turbulent 
heat transfer from various popular methods could be divergent at low values of Tw/Tt 
depending on the Mach number range. Consequently, accurate experimental data for the 
turbulent boundary layer with low ratios of wall temperature to total temperature are 
needed to determine the most acceptable prediction method for design purposes. The 
investigations of Bertram and Neal (ref. 4) and Wallace (ref. 10) give some data with low 
Tw/^t for comparison purposes, but more extensive and accurate data are necessary in 
order to establish reliability of a prediction method. 

Another perplexing problem area related to the problem of turbulent flow is the pre- 
diction of boundary-layer transition. One of the most controversial facets of high-speed 
boundary-layer transition research has been the effect of heat transfer on the transition 
Reynolds number. Various results show that heat transfer has no effect on transition 
Reynolds number (refs. 11 and 12), has a fairly strong effect on transition Reynolds num- 
ber (refs. 13 and 14), and has even caused a transition reversal (refs. 15 and 16). There 
have been numerous explanations, although without adequate substantiation, for the con- 
flicting experimental results. 

The present investigation proposes to provide additional information on two- 
dimensional, zero-pressure-gradient turbulent-boundary-layer heat transfer and transi- 
tion over a wide range of ratios of wall temperature to total temperature and to determine 
which of the theories presently considered most reliable best predicts the variation of 
surface heat transfer with wall cooling. Internally cooling a flat -plate model provided 
wall temperature ratios over the range 0.19 i T^/Tt ^ 0.7. Data were obtained for two 
unit Reynolds numbers with the stagnation temperature held constant at approximately 
533° K and for two angles of attack giving local Mach numbers of 6.0 and 4.9. In all 
cases laminar flow occurred near the leading edge, but transition always ended before 
the midpoint of the plate. 


SYMBOLS 


Cf skin friction coefficient 

Ciff specific heat of model material 


2 



M 

Nst 

P 

q 

R 

^crit 

^,Ax 

r 

T 

T» 

t 

u 

X 

Ax 

y 

a 

y 

X 


specific heat of air at constant pressure 
Mach number 

q 

Stanton number, r 

pucp(Taw - Tw) 

absolute pressure 
surface heat -transfer rate 
unit Reynolds number, ^ 
minimum critical Reynolds number 

pgUgAx 

transition-length Re 3 molds number, 

Me 

local Reynolds number based on distance from leading edge to inflection 
point (see fig. 11) 

recovery factor 

temperature 

reference temperature 

time 

velocity 

longitudinal distance from plate leading edge 
longitudinal extent of transition region 
spanwise distance from plate center line 
angle of attack 

ratio of specific heats for air 
skin thickness 



M 


viscosity 


p density 

Subscripts: 

aw adiabatic wall 

b based on distance from leading edge to beginning of transition (fig. 11) 

e local conditions at edge of boundary layer 

i local incompressible 

t stagnation 

tr based on distance from leading edge to end of transition (fig. 11) 

V based on distance from virtual origin 

w wall 

X based on distance from leading edge 

oo undisturbed free stream 

APPARATUS AND TESTS 
Tunnel 

Tests were conducted in the Langley 20-inch hypersonic tunnel. This timnel has a 
fixed nozzle block and a rectangular test section and is of the blowdown type exhausting to 
the atmosphere through a movable second minimum. An annular ejector downstream of 
the test section insures lower starting and running pressures. Normal operating stagna- 
tion pressure can be as high as 37 atmospheres, and by using electrical heaters, the stag- 
nation temperature can be as high as 589® K. Test times as long as one-half hour are 
possible. A calibration of the test core indicates that for the present test conditions, the 
Mach number is 6.02 ± 0.02. The water concentration in the tunnel airflow was measured 
to be 1 to 2 parts per million. 
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A model injection system using a pneumatic cylinder was located on top of the tun- 
nel directly above the test section. A rectangvilar opening (approximately 55.88 cm by 
36.8 cm) allowed the model to be injected to the center line of the test section after the 
tunnel flow had been established. 

Model and Instrumentation 

The flat -plate configuration used in this investigation was intended to produce a 
two-dimensional, zero-pressure-gradient boundary layer along the test surface. The 
construction and dimensions of the model are shown in figure 1 along with the orienta- 
tion of the instrumentation. The plate was constructed of AISI 405 stainless steel with the 
leading-edge thickness measured as 0.0038 ± 0.0008 cm across the span. Coolant pas- 
sages were located on each side of the model center line as shown in figure 1. For the 
present tests liquid nitrogen was injected under pressure into one passage on each side 
of the model center line and allowed to exhaust from the others. The model was mounted 
to the injection carriage with a strut attached to the model on the side opposite the instru- 
mented surface. (See fig. 1.) 

A 2.54-cm-wide cavity was milled into the back of the model until the metal skin 
between the bottom of the cavity and the flat surface was 0.076 ± 0.0025 cm thick. Thirty- 
gage iron-constantan thermocouples were then spotwelded to the undersurface of the thin 
skin along the center line of the model; see table I for the x- and y-locations of the 
thermocouples. The skin thickness was measured at each thermocouple location for use 
in the data reduction. The cavity was then filled with a foamed plastic substance with a 
very low thermal conductivity (k « 0.021 watt/meter-°K). A stainless -steel cover plate 
which faired the cavity into the bottom of the plate was then applied. 

Pressure orifices were located along a line 3.81 cm from the model center line at 
the x-locations shown in table I. The orifices were made from standard 0.23 cm outside 
diameter monel tubing, the end of which was mounted flush with the flat surface of the 
plate. The pressure and thermocouple leads were run along the base of the strut and out 
of the test section. The pressure leads were connected to 0 to 6.9 x 10^ N/m^ diaphragm- 
type electrical transducers (Statham gages). 

End plates were applied to the model in both the up and down positions as shown in 
figure 2. The leading-edge bevel angle of the end plates was 15° and the edge thicknesses 
facing the flow were approximately 0.005 cm. The end plates were designed to enclose 
the leading-edge shock at all angles of attack. 

Test Conditions and Procedures 

Most of the heat -transfer and pressure tests were made with conditions which would 
give the largest extent of turbulent-boundary-layer flow over the model. Most of the data 
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were obtained with a stagnation pressure of approximately 3.55 x 10® N/m^ at a stag- 
nation temperature of approximately 533° K corresponding to a free-stream unit 
Reynolds number of approximately 2.64 X 10® per cm. Selected heat -transfer tests 
were performed at an angle of attack of 0° with a stagnation pressure of approximately 
1.83 X 10® N/m2 and a stagnation temperature of approximately 533° K; the corre- 
sponding free-stream Reynolds number was approximately 1.46 X 10® per cm. 

For the tests at an angle of attack of 8.1° (compression), the stagnation pressure 
and temperature were 3.55 x 10® N/m^ and 533° K, respectively. These conditions 
resulted in a local unit Reynolds number of approximately 3.82 x 10® per cm and a local 
Mach number of 4.9; the local Reynolds number and Mach number were calculated by use 
of oblique shock relations, the solutions for which are tabulated in reference 17. 

The model was cooled by injecting liquid nitrogen through the coolant passages 
shown in figure 1. Since the test facility was not airtight, it was necessary to cover the 
model during the coolii^ process to prevent frost formation on the instrumented surface. 
The model, positioned in the injection housing, was covered with aluminum foil and then 
a rubber sheet connected to each side of the injection housing. When the model had been 
cooled to the desired temperature, the tunnel was started, and the model was injected 
into the established stream. As the model plunged from the cavity toward the stream, 
the rubber sheet was broken and the aluminum foil was ripped away from the model by 
the tunnel air and dissipated downstream. An observer watched the injection for each 
run to insure that no aluminum foil or surface frost was present when the model reached 
the center line of the tunnel test core. The model moved from the sheltered cavity to the 
tunnel center line in approximately 0.25 second. 

The following nominal conditions were used in calculating parameters from all 
methods; 


II 

o 

o 

Mg = M«, = 6.0 

Tt,e = Tt^co = 533° K 
Tg = Too = 65° K 

a = 8.1° 

Mg = 4.9 
Tt^g = 533° K 
Tg = 92.2° K 


The model wall temperature was variable from run to run. Typical surface temperature 
distributions at t = 0 are presented in figure 3. The wall temperature generally varies 
less than ±15° about a mean value; this mean value has been used for the calculations. 
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DATA REDUCTION 


All the experimental data were reduced by use of a nominal free-stream Mach 
number of 6.0 and measured values of free-stream stagnation pressure and temperature. 

Pressure Data 

The electrical outputs from the pressure transducers were recorded on a high- 
speed digital readout recorder about 5 seconds after the model was in final position in 
the test section of the tunnel. From the calibration data for each transducer, electrical 
outputs were converted to pressure readings on a card-programed computer. The accu- 
racy of the transducers was 0.25 percent of full-scale pressure. 

Heat-Transfer Data 

The electrical outputs from the thermocouples were recorded on a high-speed 
digital readout recorder. The signal from each thermocouple was sampled 20 times 
each second, converted to a binary digital system, and recorded on magnetic tape. When 
the model was first positioned at the center line of the tunnel test section, 1 second of the 
temperature -time data was fitted by the method of least squares to a second-degree poly- 
nomial of the form 

T^ = a + bt + ct^ (1) 

where a, b, and c are constants. The time derivative of temperature used to calcu- 
late the heating coefficients was computed at the first point of the curve fit. (The deriva- 
tive was constant for at least the first five points.) 

The model wall temperature at the start of the tests varied from approximately 
100° K to 341° K, Because of the quick insertion into the test flow (approximately 
0.25 sec), the model was considered to have been subjected to a step function in aerody- 
namic convective heat input. In the absence of radiative and conductive heating, the local 
surface heating rate for the model was expressed as 

4 = (2) 

where the model density p^ was 7700 kg/m^. The specific heat of the model material 
is shown in figure 4 as a function of temperature. These measurements of 
were made by the National Bureau of Standards with a test specimen from the present 
model and are accurate to within 0.3 percent for T^ less than 100° K and 0.1 percent 
for T^ greater than 100° K. The skin thickness X was measured at each thermo- 
couple location and was 0.076 ± 0.0025 cm. 
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Radiative heating for the most extreme conditions encountered in this investigation 
(tunnel wall temperature ~ 333° K, model wall temperature « 100® K) was calculated to 
be less than 0.50 percent of the aerodynamic heating and thus was neglected for all data. 
Conduction heating was calculated by use of the three-point finite-difference method 
described in reference 18. Because spanAvise and chordwise conduction heating was cal- 
culated to be generally less than 1 percent of the convection heating, no corrections for 
conduction were applied to the experimental data. The heat conducted to the insulatii^ 
material backing the instrumented surface was calculated to be less than 1 percent of the 
convective heat input and was therefore neglected. 

The Stanton number was calculated from the relations: 

^Stjoo “ T T! ~ TTjj; _ rp Y 

’ PooUooCp(^iaw ■‘^wj 
or 

q . 

Pe^ie^pC^aw “ 

where 

Taw = Te(l+r2i^Me2j (5) 

The recovery factor r was taken to be 0.845 for laminar flow and 0.89 for turbulent 
flow. 


RESULTS AND DISCUSSION 

A schlieren photograph of the model without end plates alined at o; = 0® (Moo = 6) 
to the free stream is shown in figure 5. The outline of the model and rectangular strut 
is indicated on the figure. The apparent thickness of the leading-edge shock wave is due 
to a small roll angle of the model as well as some diffraction. The reflection of the 
leading-edge shock wave from the tunnel-wall boundary passes the plane of the model 
well downstream of the rear of the model and thus causes no interference on the instru- 
mented surface. 

Selected measurements showing the effects of wall cooling on the local surface pres- 
sure distribution were obtained and are shown in figure 6. The scatter in the data occurs 
because the pressures had not completely "settled out" when the data were recorded. 
Pressures were recorded approximately 5 seconds after model injection so that the model 
surface would be at a low and uniform temperature. If the pressures had been allowed 
to settle, the accuracy would have been better but large temperature gradients would have 
existed along the plate. When this procedure was used, the pressures were as much as 
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6 percent higher than the inviscid prediction and there was no trend in the data scatter. 

By comparing part (a) with part (b) and part (c) with part (d) of figure 6, the data indicate 
that wall cooling has no obvious effect on the pressure distribution at either angle of 
attack; however, wall temperature, in the absence of viscous induced effects, should have 
no effect on the pressure distribution and thus local flow conditions. For all calculations 
of local flow conditions, the local static pressure was assumed to be the inviscid value. 

Extraneous Effects on Heat Transfer 

Spanwise and end plate effects .- At two locations on the plate (x = 29.84 cm and 
36.50 cm), two spanwise thermocouples were located 0.635 cm on each side of a thermo- 
couple on the instrumented center line. For all tests, the end of transition along the 
center line had occurred before x = 29 cm; therefore, the spanwise instrumentation was 
always located in the turbulent regime. The calculated heating rates for the three span- 
wise thermocouples at a given x-station generally agreed within 5 percent at an angle of 
attack of 0° and within 4 percent at an angle of attack of 8.1®. The variation of the heating 
rates across the short span (1.27 cm) was random, was not affected by wall temperature, 
and was within the accuracy of the data. Specific data points for the spanwise instrumen- 
tation are omitted from all figures in the interest of clarity. In general, the values of 
spanwise heating for turbulent flow would fall within the bounds of the symbol for the 
center-line heating at the same x-location. 

Effects of end plates on the surface heating (Stanton number) are shown in figure 7 
at aisles of attack of 0® and 8.1® for a ratio of wall temperature to total temperature of 
0.60. A comparison of the heating data for no end plates, end plates up, and end plates 
down (see fig. 2 for orientation) indicates that the addition of end plates has no appreci- 
able effect on the surface heating along the length of the plate at either angle of attack. 

The width of the plate was such that disturbances moving from the corner of the leading 
edge toward the center line of the plate along a Mach line originating at the leading edge 
would not reach the center line. The data confirm that the heat -transfer distribution 
down the center portion of the plate was independent of the end plate orientation. The 
heating data shoAwn in subsequent figures were obtained with the "end plates down" con- 
figuration. End plates were used only as a safeguard against extraneous disturbances 
feeding onto the plate and were used in the down position because the instrumented sur- 
face could easily be covered during the coolii^ process. 

Frost formation effects .- The instrumented surface of the plate as well as the 
leading edge of the model was shielded while the plate was cooled. When the plate was 
injected, visual checks were made for all tests to insure a frost -free surface during the 
data-acquisition period. Of approximately 80 runs made during this investigation, slight 
frost appeared for approximately 20 percent of the runs. Repeat nms, one with slight 
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surface frost and one with no frost, showed that slight frost had little effect on the sur- 
face heating and transition location. An example of the effect of heavy frost formation 
on the surface heating and transition location is shown in figure 8. Heavy frost on the 
model surface was present when the plate was cooled while exposed to the atmosphere. 
The wall temperature distributions for the heavy-frost and no-frost cases are shown to 
be similar in the upper part of the figure. The heat -transfer distribution with heavy 
frost (shown by the shaded symbols) shows much data scatter and is at a lower level than 
the no-frost heating distribution. Heavy frost delays boundary-layer transition as evi- 
denced by the location of peak heating. This increased length to the end of transition 
might be considered unusual since many investigators have shown that surface roughness 
usually promotes transition; however, as shown in reference 19, small roughness may 
delay transition. The sublimation of frost at the model surface may also affect the tran- 
sition location. In light of the observed effects of frost formation, only data for which 
the instrumented surface was frost free are presented hereafter. 

Effects of Reynolds Number and Wall Cooling on Surface Heating 

Tsrpical examples of the surface heat -transfer distributions are shown in figure 9 
for an angle of attack of 0° and in figure 10 for an angle of attack of 8.1°. The free- 
stream Stanton number is plotted against the distance from the leading edge of the plate 
for various ratios of wall temperature to total temperature. Several runs are plotted 
for each ratio of wall temperature to total temperature to show the repeatability of the 
heat -transfer distributions and transition location. Heat -transfer data are repeatable 
within 10 percent for all cases, and the transition location is repeatable to within 2 cm. 
Heat -transfer data for representative ratios of wall temperature to total temperature 
are provided in tables H and HI. 

At an angle of attack of 0° (Re/cm = 2.64 x 10®), the laminar heat -transfer data 
prior to transition are well predicted by the Monaghan T' method (ref. 20). A compari- 
son of parts (a) , (b) , and (c) of figure 9 shows that the laminar heating is a weak function 
of the ratio of wall temperature to total temperature for 0.2 ^Tw/Tt^O .6 as predicted 
by the T' method. No laminar data are presented for an angle of attack of 8.1° (fig. 10) 
since transition occurred before the first thermocouple in every case. Five runs were 
made at an ai^le of attack of 0° and Re/cm = 1.46 X 10^ which are not presented in 
figure 9. The laminar data for these rims were also predicted well by Monaghan's 
T' method. 

The virtual origin of the turbulent boundary layer was chosen as the location of 
peak heating near the end of transition (as in ref. 4) and the Reynolds number based on 
this origin (Re,v) illustrated in figure 11. The variation of the turbulent heat transfer 
with turbulent Reynolds number (Re,v) based on the distance from the virtual origin is 
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presented in figure 12 for Mg = 6.0 and in figure 13 for Me = 4.9. The heating data 
are presented as the ratio Nst^e/^Stji where the incompressible Stanton number^ was 
obtained from reference 21 at the same Reynolds number as N 3 ^^g. Generally, several 
runs are presented for each ratio of wall temperature to total temperature including some 
data for Rg/cm = 1.46 x 10®. At the lowest ratio of wall temperature to total tempera- 
ture for both Mg = 6.0 and 4.9, the experimental Stanton number decreases more slowly 
with Reynolds number than the incompressible Stanton number (shown by the increase of 
Nst,e/^St,i with Re,v)* As the wall temperature is increased, Ng^^gy^N 3 ^^j becomes 
constant with Re,v until at the higher wall temperatures Ng^^e/^St,! decreases weakly 
with Reynolds number. Therefore, the slope of the heat transfer with Reynolds number 
is a function of the ratio of wall temperature to total temperature. Note that the data for 
Re^v l®ss than approximately 10® generally fall below the data for Re,v greater than 
106 . A similar effect was found in reference 4 for a large body of turbulent heating data. 

Also presented in figures 12 and 13 are the predictions for the turbulent heating 
from three methods: the Monaghan turbulent T’ method (ref. 5), the Spalding and Chi 
method (ref. 2), and the method (ref. 9). The skin-friction predictions from the 

T’ (using the Karm^ and Schoenherr skin-friction equation) and Spalding and Chi methods 
were converted to heat transfer by using a modification of Karma’s Reynolds analogy as 
presented in reference 4. At Mg = 6 (fig. 12), the method which best predicts the level 
and trend of the data through the range of ratios of wall temperature to total temperature 
is the Spalding and Chi method. The method consistently underpredicts the level 

of the data for each ratio of wall temperature to total temperature and does not match the 
trend of the data with increasing Rg y as well as the Spalding and Chi method. The 
T' method substantially overpredicts the level of the data at the lower ratios of wall tem- 
perature to total temperature; as the ratio of wall temperature to total temperature is 
increased to 0.6 and above, good agreement exists between the T' prediction and the 
experimental data. As previously mentioned, comparisons were with data for which 
Re,v was greater than 10®. At Me = 4.9 (fig. 13), the data fall between the predic- 
tions of the Spalding and Chi and Pj.Pr methods but generally favor the level of the 
Spalding and Chi method. The T' method again overpredicts the heat transfer at low 
ratios of wall temperature to total temperature. 

To illustrate clearly the effect of wall cooling on turbulent heat transfer, the data 
are presented as Stanton number ratio as a function of wall cooling for a specific value 
of Re,v in figure 14 for Me = 6 and in figure 15 for Me = 4.9. Each data point in 
figures 14 and 15 was obtained from a fairing of Ng^. e/Nst,i against Re,v at the 

2 Ngj^j^ was obtained by multiplyii^ the skin-friction coefficient from the Karm^- 

Schoenherr imcompressible equation by Karman's value of Reynolds analogy at the same 
Reynolds number as Ng^^g. 
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particular value of Re v indicated in figure 14 or 15. In this form the data show that 
at both Mach numbers the turbulent heat transfer is approximately independent of the 
wall cooling for the range of ratios of wall temperature to total temperature for this 
investigation. Predictions by the same methods shown in figures 12 and 13 as well as 
predictions by the Van Driest method (ref. 22) using a Reynolds analogy factor of 1.16, 
the Sommer and Short T’ method (ref. 7), the Eckert T’ method (ref. 6), and the Winkler 
and Cha method (ref. 8) are included in figures 14 and 15 for comparison with the data. 
These engineering methods are most commonly used to estimate hypersonic turbulent 
boundary-layer skin friction and heat transfer. At Me = 6, the prediction from the 
Spalding and Chi method agrees with both the level and the trend of the data with wall 
cooling; at Mach 4.9 the Spalding and Chi method predicts the trend with wall cooling but 
slightly overpredicts the level of the turbulent heating. The Pj-Mr method predicts the 
trend but underpredicts the magnitude of the data at both Mach 6 and 4.9. Clearly, at 
both Mach numbers, the Van Driest and T' methods substantially overpredict the heat 
transfer at the low ratios of wall temperature to total temperature, and the Winkler and 
Cha method substantially underpredicts the data at the same low ratios of wall tempera- 
ture to total temperature. 

Previous investigations by Bertram and Neal (ref. 4), Wallace (ref. 10), and Nerem 
and Hopkins (ref. 23) have shown that the Spalding and Chi method best predicts their tur- 
bulent skin-friction and heat -transfer results for the lower ratios of wall temperature 
to total temperature (or enthalpy). The results of Wallace (Me = 4 to 9) and Nerem and 
Hopkins (Me = 2.5 to 3.5) were obtained in shock tunnels where low ratios of wall enthalpy 
to total enthalpy were achieved on a room-temperature model by elevating the stagnation 
enthalpy of the stream; in contrast, the model was cooled and the stagnation temperature 
held constant for the results of this investigation. The resultii^ comparisons of experi- 
mental data with turbulent prediction methods yielded the equivalent result that the 
Spalding and Chi method best predicted the data. 

A recent investigation by Hopkins, et al. (ref. 24) indicated that the methods of 
Coles (ref. 25) and Van Driest (ref. 22) best predict the turbulent skin-friction data avail- 
able above Mach 4 or 5 and Tw/T^w ^ when the predictions are based on a momen- 
tum thickness Reynolds number rather than on a length Reynolds number. Turbulent 
heat -transfer data were also better predicted by these same methods when a Re 3 molds 
analogy factor ^2Nst/Cf^ of 1.0 was used instead of the extension of Von Karman’s 
Reynolds analogy as used herein and when the predictions were based on an energy thick- 
ness Reynolds number (obtained by integrating the surface heat -transfer distribution) 
instead of a length Reynolds number. The present turbulent heating data, included in the 
analysis of reference 24, were well predicted by Coles' and Van Driest' s methods and 
were significantly underpredicted by Spalding and Chi's method when an energy thickness 
Reynolds number and a Reynolds analogy factor of 1.0 were used. 
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In reference 24 a Reynolds analogy factor of 1 was experimentally measured for 
Mach 6.8 and 7.4 ® and was assumed to remain equal to 1.0 at other Mach 

numbers and ratios of wall temperature to adiabatic wall temperature for comparing data 
with prediction methods. If a Reynolds analogy factor of 1.0 had been used with the pre- 
diction methods presented in this report instead of the modified Karman Reynolds analogy, 
the predictions would be reduced by 9 to 17 percent. A compilation of available experi- 
mental measurements of Reynolds analogy for compressible, turbulent -boundary- layer 
flow presented in reference 26 tend to favor Karma's values, in general, for T^/Tt 
greater than 0.4. Predictions are also sensitive to the choice of the virtual origin for 
the turbulent boundary layer. For the data presented in figure 9, the momentum thick- 
ness Reynolds number was obtained as in reference 24 and this momentum thickness 
Reynolds number extrapolated to zero by using the Spalding and Chi method. By using 
this new virtual origin (zero-momentum-thickness location), the Spalding and Chi predic- 
tions for a given x-station (x > 28 cm) decreased from 7 to 18 percent depending on the 
Resmolds number and the ratio of wall temperature to total temperature. It thus appears 
that selecting a "best" method for predicting hypersonic turbulent heating will depend 
upon the methodology and assumptions involved in the application of the methods (choice 
of Reynolds analogy factor, virtual origin, and Reynolds number). Data and theory com- 
parisons at very high Reynolds numbers, such as represented by Wallace’s results (see 
presentation of his data in ref. 27), are little affected by at least the choice of virtual 
origin. 


Effects of Wall Cooling on Boundary-Layer Transition 

Boundary-layer transition Reynolds numbers were determined from the heat- 
transfer distributions as indicated in figure 11. The Reynolds number for the beginning 
of transition was taken as the local unit Reynolds number times the x-distance from the 
leading edge of the plate to the minimum Stanton number; the Reynolds number for the 
end of transition was taken as the local unit Reynolds number times the x-distance from 
the plate leading edge to the location of the maximum Stanton number. (See fig. 11.) 

The effect of wall cooling on the transition Reynolds number for the beginning and 
end of transition at Mach 6 and 4.9 are shown in figure 16 and these data are given in 
table IV. At Mach 6 transition results are presented for Re/cm = 2.65 x 10® along with 
limited results for Re/cm = 1.46 x 10®. The Mach 6 results for Re/cm = 2.65 x 10® 
show that the transition Reynolds number increases approximately 40 percent for a 
decrease in the ratio of wall temperatvire to adiabatic wall temperature from 0.7 to 0.2. 
For Rgyfcm = 1.46 X 10^, the transition Re5molds number increases approximately 60 per- 
cent for a Ty^Tg^^ decrease from 0.7 to 0.2. The Reynolds number for the beginning 
and end of transition increase approximately the same percentage through the range of the 
ratio of wall temperature to adiabatic wall temperature at a given unit Reynolds number. 
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At Mach 4.9, the Reynolds number for the end of transition increased approximately 
30 percent for a '^w/'^aw decrease from 0.7 to 0.2. The Reynolds number for the 
beginning of transition at Mach 4.9 is not shown because transition had begun before the 
first measuring station for each case. 

A comparison of fairings of the present transition data with faired results from 
previous investigations using sharp-leading-edge (leading-edge thickness ^ 0.0051 cm) 
flat plates (refs. 11, 15, 28, 29, and 30) is shown in figure 17 for Mach numbers from 2.4 
to 10.2. There appears to be no consistent variation of transition Re 3 molds number with 
decreasing ratio of wall temperature to adiabatic wall temperature. Depending upon the 
investigation, the transition Reynolds number increases, decreases, or first increases 
and then decreases with decreasing Tw/'^aw* The reversal trend for the Mach 8.2 data 
of Richards and Stollery (ref. 15) was characterized by the transition location moving 
rearward on the plate with wall cooling until the boundary layer was laminar over the 
length of the plate; further cooling caused transition to reoccur on the plate and subse- 
quently to move rearward again beyond the end of the plate. This transition-reversal 
phenomenon has also been observed to occur on cones at lower ratios of wall temperature 
to adiabatic wall temperature. (See ref. 16.) 

Since there is no theoretical approach to explain the effects of variations of flow 
variables and boundary conditions on compressible boundary-layer transition, experi- 
mentally observed transition results have not been satisfactorily correlated. Some 
insight into the reasons for the differing experimental results may be available from a 
cursory look at linearized stability theory if one could assume that the variation of the 
minimum critical Reynolds number (that Reynolds number below which disturbances in 
the boundary layer will not be amplified) is generally similar to the variation for the 
transition Re}molds number. This correspondence of Rcrit transition Resmolds 

number may not exist. Mack's recent solution to the compressible stability equations 
(ref. 31) has indicated that numerous modes of amplified solutions to the viscous stability 
equations occur. Solutions prior to Mack's solution were concerned only with the first 
mode. The numerical results show that wall cooling is stabilizing to the first mode but 
destabilizes the second and higher modes; also wall cooling increases the value of the 
frequency which is most unstable in the boundary layer. Thus it appears that depending 
on the spectrum of disturbances in the boundary layer as well as the Mach number and 
Reynolds number, wall cooling may or may not produce a stabilizing effect on transition. 

If the variation of transition Rejmolds number with wall cooling does in fact depend 
upon the disturbances present in the flow field, it is not surprising that the trends of the 
experimental results with wall cooling in figure 17 do not agree; each of the sets of data 
were obtained in different facilities which, in general, have a different disturbance spec- 
tnim in the free stream. The complex interrelationships between the various modes of 
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amplification of disturbances as well as the dependency of boundary -layer stability upon 
the frequencies of disturbances present combine to form a structure which could in prin- 
ciple allow any variation of transition Reynolds number with wall cooling. The further 
implication of the theoretical results is that since the flow-field disturbances found in 
flight would, in general, be different from those found in wind tunnels, the effect of wall 
cooling on the transition Reynolds number found in wind-tunnel tests may not be repre- 
sentative of the effect found in flight tests. It is highly probable that in order to obtain 
a meaningful correlation of wind-tunnel transition results, additional parameters such as 
the spectrum of free-stream disturbances will be needed. 

The variation of transition Reynolds number with unit Reynolds number for selected 
values of the ratio of wall temperature to total temperature is shown for Mach 6 on the 
left-hand side of figure 18. The transition Reynolds number increases with the local unit 
Reynolds number to the 0.4 to 0.45 power through the wall temperature range. In con- 
trast, the data of Richards and Stollery (ref. 15) at Mach 8.2 shown on the right-hand side 
of figure 18 show that the transition Reynolds number increases much more rapidly with 
unit Reynolds number as the ratio of wall temperature to total temperature is decreased. 
There are many references available which analyze the effect of unit Reynolds number on 
the transition Reynolds number for supersonic and h 3 T»ersonic speeds (see, for example, 
ref. 32), but these analyses are only for moderate values of the ratio of wall temperature 
to total temperature. Generally, previous investigators found the unit Reynolds number 
effect to be similar to the results for the present investigation. 

Information concerning the effect of wall cooling on the extent of the transition 
region is available from the data at Mach 6. Potter and Whitfield (ref. 32) have shown 
for adiabatic wall and near -adiabatic wall conditions that the Reynolds number based on 
length of the transition region for two-dimensional flow is fairly independent of unit 
Reynolds number and leading-edge geometry. The correlation which resulted was of the 
form 

^,Ax ~ ^(^e>^,infl) (®) 

where these Reynolds numbers are illustrated in figure 11. This correlation was shown 
to be independent of sweep angle at Mach 8 in reference 11. The present data are pre- 
sented according to the Potter and Whitfield correlation in figure 19. The solid lines 
shown in figure 19 for each Mach number are taken from the Potter and Whitfield correla- 
tion. Data from the present investigation representing a variation in the ratio of wall 
temperature to total temperature from 0.19 to 0.70 are in substantial agreement with the 
correlated curves. The agreement indicates that at least for Mach 6, the Potter and 
Whitfield correlation for the length of the transition region is independent of heat transfer 
as well as unit Re 3 molds number and leading-edge geometry. 
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CONCLUSIONS 


A study of the effect of wall cooling on zero-pressure-gradient turbulent-boundary- 
layer heat transfer and laminar-boundary-layer transition has been conducted at a free- 
stream Mach number of 6 and Reynolds numbers per cm of 2.64 X 10® and 1.46 X 10® in 
the Langley 20 -inch hypersonic tunnel. The wall temperature of a sharp-leading-edge 
flat plate was varied by internal cooling to yield a range of the ratio of wall temperature 
to total temperature from 0.19 to 0.70 (with no surface frost) for a total tempera- 

ture of 533° K. With the model positioned in the tunnel test section at angles of attack of 
0° and 8.1°, the local Mach numbers were 6.0 and 4.9, respectively. The conclusions 
derived from this study are as follows: 

(1) Although light frost had little effect on the results, heavy frost formation on the 
model significantly affected the level of the surface heating as well as the location of 
boundary-layer transition; therefore only data for frost -free conditions were included in 
the analysis. 

(2) Decreasing the ratio of wall temperature to total temperature (o.l9 ^ Tw/Tt £ 0.7) 
had little effect on the heat -transfer coefficient at either Mach 6 or Mach 4.9. 

(3) With the virtual origin chosen as the point of peak heating near the end of transi- 
tion, the Spalding and Chi method gave a good prediction of the turbulent heat transfer (for 
a local Reynolds number based on the distance from the virtual origin greater than 106) 

at both Mach 6 and Mach 4.9 over the entire range of the ratio of wall temperature to total 
temperature. The turbulent T' and Winkler and Cha methods were found to be inadequate 
for predicting the heating for the low ratios of wall temperature to total temperature. 

This conclusion is true only when it is assumed that the virtual origin of turbulent flow is 
at peak heating and that the modified K§.rm&n Reynolds analogy applies. 

(4) For a given Mach number and unit Reynolds number, transition Reynolds number 
increased 30 to 60 percent with decreasing ratio of wall temperature to total temperature, 
and no transition reversal occurred. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 3, 1970. 
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TABLE I.- LOCATIONS OF INSTRUMENTATION 
[Dimensions are in cn^ 


Thermocouple 

X 

y 

Thermocouple 

X 

y 

1 

7.30 

0 

20 

29.84 

-0.635 

2 

9.21 



21 

30.80 

0 

3 

10.16 



22 

31.74 



4 

11.12 



23 

34.60 



5 

13.03 



24 

35.57 



6 

13.97 



25 

36.50 



7 

14.92 



26 

37.48 



8 

15.88 



27 

37.48 

.635 

9 

16.84 



28 

37.48 


635 

10 

19.38 



29 

38.41 

0 

11 

21.27 



30 

39.38 



12 

22.22 



31 

42.22 



13 

23.19 



32 

43.19 



14 

24.12 



33 

44.14 



15 

26.99 



34 

45.10 



16 

27.94 



35 

46.04 



17 

28.89 



36 

47.00 



18 

29.84 


37 

49.20 



19 

29.84 

0.635 

38 

50.18 




Pressure 

orifice 

X 

1 

7.62 

2 

12.70 

3 

17.78 

4 

35.57 

5 

40.00 

6 

45.73 

7 

50.80 


y 

3.81 
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TABLE n.- HEAT-TRANSFER DATA FOR Mp = 6.0 


[m^ = 6.0; Q! = 0^ 


Re/cm = 0.31 X 10°; 


Nst,oo 

Tw/Tt = 0.7 
2.68 X 10^ 


Re/cm = 0.26 x 10®; 


10.16 

11.11 

13.02 

13.97 

14.92 

15.87 

16.84 

19.38 
21.27 
22.22 

23.19 

24.13 
27.94 
28.90 

29.85 
31.75 
34.61 
35.58 
36.51 
37.49 
38.41 

39.38 
42.21 

43.19 

44.13 
45.10 
46.04 
47.00 

49.20 


TABLE II.- HEAT -TRANSFER DATA FOR Me = 6.0 - Continued 


X, cm 

^StjOO 

X, cm 


Re/cm = 0.26 X 10®; T^/Tt = 0.5 

Re/cm =s 0.27 X 10^; T^/Tt = 0.4 

7.30 

2.98 X 10^ 

7.30 

2.95 X 10^ 

9.20 

3.12 

9.20 

2.98 

10.16 

3.25 

10.16 

2.83 

11.11 

3.91 

11.11 

3.03 

13.02 

4.65 

13.02 

3.13 

13.97 

6.56 

13.97 

4.67 

14.92 

7.69 

14.92 

6.01 

15.87 

8.69 

15.87 

7.19 

16.84 

8.54 

16.84 

7.65 

19.38 

8.54 

19.38 

8.54 

21.27 

7.95 

21.27 

7.71 

22.22 

7.92 

22.22 

7.72 

23.19 

7.51 

23.19 

7.57 

24.13 

7.21 

24.13 

6.86 

27.94 

6.67 

27.94 

6.77 

28.90 

6.82 

28.90 

6.26 

29.85 

6.75 

29.85 

6.59 

31.75 

6.95 

31.75 

6.25 

34.61 

6.46 

34.61 

6.13 

35.58 

6.30 

35.58 

6.40 

36.51 

6.16 

36.51 

5.98 

37.49 

6.02 

37.49 

5.63 

38.41 

6.14 

38.41 

5.64 

39.38 

5.89 

39.38 

5.17 

42.41 

5.94 

42.21 

5.63 

43.19 

5.58 

43.19 

5.35 

44.13 

5.57 

44.13 

5.43 

45.10 

5.64 

45.10 

5.25 

46.04 

5.46 

46.04 

5.02 

47.00 

5.42 

47.00 

4.94 

49.20 

5.40 

49.20 

4.79 
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TABLE II.- HEAT-TRANSFER DATA FOR Mg = 6.0 - Continued 




TABLE n.- HEAT-TRANSFER DATA FOR Mg = 6.0 - Continued 


X, cm 

8 

X, cm 

Nst,oo 

Re/cm = 0.14 X lO^; Tw/Tt = 0.6 

Re/cm = 0.14 X 10®; T^/Tt = 0.4 

7.30 

3.65 X 104 

7.30 

3.65 X 1q4 

9.20 

3.69 

9.20 

3.20 

10.16 

3.62 

10.16 

2.78 

11.11 

3.53 

11.11 

3.07 

13.02 

3.91 

13.02 

3.42 

13.97 

4.28 

13.97 

3.95 

14.92 

4.73 

14.92 

4.30 

15.87 

5.45 

15.87 

4.47 

16.84 

5.83 

16.84 

4.63 

19.38 

7.81 

19.38 

5.77 

21.27 

9.04 

21.27 

7.09 

22.22 

9.66 

22,22 

8.10 

23.19 

9.23 

23.19 

8.45 

24.13 

8.91 

24.13 

8.51 

27.94 

8.02 

27.94 

8.42 

28.90 

7.95 

28.90 

8.37 

29.85 

7.78 

29.85 

8.27 

31.75 

8.03 

31.75 

7.82 

34.61 

7.40 

34.61 

7.23 

35.58 

7.88 

35.58 

7.92 

36.51 

6.84 

36.51 

6.86 

37.49 

6.72 

37.49 

7.05 

38.41 

6.86 

38.41 

7.27 

39.38 

6.64 

39.38 

6.53 

42.21 

6.92 

42.21 

6.63 

43.19 

6.74 

43.19 

6.49 

44.13 

6.12 

44.13 

6.47 

45.10 

6.22 

45.10 

6.56 

46.04 

5.99 

46.04 

6.16 

47.00 

6.13 

47.00 

6.03 

49.20 

6.62 

49.20 

6.15 
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TABLE m.- HEAT -TRANSFER DATA FOR Mg = 4.9 
[Mco = 6.0; a = 8.1^ 


X, cm 



X, cm 

Nst,oo 

Rg/cm = 0.39 X 10®; T^/Tt = 0.6 

Rg/cm = 0.39 X 10®; Tw/Tt = 0 

7.30 

10.82 X 10^ 


7.30 

8.93 X 10^ 

9.20 

17.67 


9.20 

13.74 

10.16 

19.55 


10.16 

16.13 

11.11 

18.18 


11.11 

17.25 

13.02 

17.41 


13.02 

17.43 

13.97 

17.23 


13.97 

17.64 

14.92 

16.63 


14.92 

17.11 

■ 15.87 

16.38 


15.87 

16.80 

16.84 

14.82 


16.84 

14.91 

19.38 

14.57 


19.38 

15.20 

21.27 

14.29 


21.27 

14.58 

22.22 

14.70 


22.22 

14.95 

23.19 

14.04 


23.19 

14.60 

24.13 

13.47 


24.13 

13.85 

27.94 

12.93 


27.94 

13.42 

28.90 

12.85 


28.90 

13.57 

29.85 

13.10 


29.85 

13.77 

31.75 

12.97 


31.75 

13.82 

34.61 

12.70 


34.61 

13.36 

35.58 

13.03 


35.58 

13.64 

36.51 

12.15 


36.51 

12.88 

37.49 

12.34 


37.49 

12.76 

38.41 

12.11 


38.41 

12.79 

39.38 

11.95 


39.38 

12.33 

42.21 

11.91 


42.21 

12.32 

43.19 

11.73 


43.19 

11.98 

44.13 

11.65 


44.13 

11.89 

45.10 

11.72 


45.10 

12.10 

46.04 

11.36 


46.04 

11.84 

47.00 

10.98 


47.00 

11.76 

49.20 

11.59 


49.20 

^ 12.08 







TABLE IV.- TRANSITION DATA 


12.70 

11.68 

12.20 

12.70 

12.20 

10.92 

9.90 
10.41 
10.92 
10.16 
"0.41 

1.17 

0.67 

0.41 

9.40 

9.91 
9.14 
0.15 
9.40 


xj-j., cm 

21.59 

21.09 

20.32 

18.80 

18.80 

16.77 

17.78 

16.36 
16.51 

17.12 
17.15 
15.89 
14.74 
14.74 

16.13 

15.37 
14.49 
16.51 
14.61 


Tw/T 

t Re/cm 

xtj, cm 

Xj-J., cm 


Me = ( 

3.0 



0.580 

.600 

.700 

.190 

.400 

.550 

.580 

.600 


2.64 X 10^ 

2.70 

3.08 

1.39 

1.44 

1.46 


9.90 

9.90 

8.13 

18.04 
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Figure 1.- Model dimensions and instrumentation. (All dimensions are in cm.) 









Figure 5.- Schlieren showing flow details and model orientation. Me = 6.0; Re/cm = 2.64 x 10^. L-70-1622 
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Figure 7.- Effect of end plate orientation on the surface heat transfer. Tw/Tf = 0.60. 
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Figure 9.- Heat-transfer distributions for various ratios of waii temperature to total temperature for 
Me = 6.0; a = OO; Rg/cm = 2.64 x 105. 


37 



Figure 10.- Heat-transfer distributions for various 


3- runs 
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Figure 17.- Experimental transition results for flat plates with wall cooling. 
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Figure 18.- Effect of unit Reynolds number on the transition Reynolds number for various wall temperature levels. 
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